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The combination of the local-density approximation 
(LDA) with the rotationally invariant slave-boson the- 
ory (RISE) is used to investigate the realistic corre- 
lated electronic structure of Sr3Ru207. From Wannier- 
downfolding the low-energy band structure to a three- 
band model for the Ru(t2g) states, the interacting prob- 
lem is solved including intra- and inter-orbital Hubbard 
terms as well as spin-flip and pair-hopping interactions. 
Therewith it is possible to obtain valuable insight into the 
orbital occupations, relevant local spin multiplets and the 
fermiology with increasing correlation strength. 



1 Introduction Strongly correlated electron systems 
are not only a fascinating research field from the viewpoint 
of fundamental research, but also become more and more 
of technological relevance in the context of the demand 
for specific high-responsive behavior The very subtle bal- 
ance between the kinetic -energy term and the Coulomb re- 
pulsion among the electrons in these systems indeed of- 
ten result in e.g. absorbing magnetic, superconducting or 
thermoelectric properties. Concerning the theoretical de- 
scription on an atomistic level, the standard Kohn-Sham 
(KS) band-theory representation of density functional the- 
ory (DFT) for the solid state is usually inappropriate for re- 
alistic materials with strong electronic correlations. For in- 
stance, the Mott-insulating state as an electron localization 
in real space amounts to a complete breakdown of the con- 
ventional band-structure concept, e.g., given by the local- 
density approximation (LDA) to DFT. However the sim- 
ple LDAh-U extension for the strongly correlated regime 



Besides generic correlation-induced band-narrowing and 
-shifting, an intriguing quasiparticle structure close to 
the Fermi level is found in the neigborhood of the no- 
torious 72 pocket in the Brillouin zone. Along the F—X 
direction in k-space, that structure appears very sensi- 
tive to electronic self-energy effects. The subtle sensitiv- 
ity, connected also its manifest multi-orbital character, 
may put this very low-energy structure in context with 
the puzzling metamagnetic properties of the compound. 
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is generally not well-suited for metals, since the correlated 
low-energy behavior close to the Fermi level ep cannot be 
captured within this static mean-field method. 

Here we want to show that the combination of LDA 
with the slave-boson theory |[Tl|2l in its rotationally- 
invariant form ||3]|4] (RISE) at saddle-point may be an 
efficient approach to model strongly correlated metals. 
This LDAh-RISE scheme can account for many details of 
the quasiparticle (QP) physics close to ep of strongly inter- 
acting multi-orbital systems with modest effort compared 
to, e.g., numerically expensive Quantum-Monte-Carlo im- 
purity computations within the framework of dynamical 
mean-field theory (DMFT). 

In this work the electronic correlations effects in the 
complex Sr3Ru207 compound shall be discussed within 
the LDAh-RISE method. The puzzling ruthenate has gath- 
ered much attention in recent years due to its metamagnetic 
behavior at low temperatures, which may be related to a 
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quantum critical point ||5]|6]|2l- The paramagnetic Fermi- 
liquid displays one of the highest electronic specific heat 
values among oxides and has definitely to be placed in 
the strongly correlated regime. Because of its very del- 
icate low-energy physics, the slave-boson technique, fo- 
cussing on the strongly renormalized QP part, is an ade- 
quate framework for investigations beyond standard LDA 
approaches. 

2 The Sr3Ru207 compound The n=2 case of 
the multilayer Ruddlesden-Popper strontium ruthenates 
Sr„+iRu„02ri+3 serves as an interesting application of 
our theoretical methodology. The perovskite end member 
SrRuOs (?i=cx}) of this family is at low temperature T 
a ferromagnetic metal, whereas the also itinerant single- 
layered Sr2Ru04 compound is paramagnetic at ambient 
T and becomes superconducting below Tc^4.2 K, with 
a widely believed triplet pairing |[8] . The crystal structure 
of Sr3Ru207 in with lattice parameters a=6=5.5006A 
and c=20.725A is based on the orthorhombic space group 
Bbcb (#68) (see Fig. [U and does not show the fourfold 
symmetry of the simpler Sr2Ru04 structure. It consists of 
Ru02 bilayers, whereby the RuOg octahedra are rotated 
by 6.8°. The unit cell exhibits four Ru ions, all equivalent 
by symmetry. 

Many experiments show that the Sr3Ru207 compound 
is located just at the paramagnetic-to-ferromagnetic tran- 
sition boundary. It is still paramagnetic in zero magnetic 
field, but positioned very close to the ferromagnetic in- 
stability ifTOl . The system shows strongly renormalized 
Fermi-liquid behavior with however a strong resistivity 
anisotropy Pc/ Pab^^OQ ifTOl . With applied magnetic field 
H, the systems shows metamagnetic behavior, i.e., a very 
large dM/dH, around i/=5.5(7.7)T for H\\abic) fUl . 
Furthermore this metamagnetic region may be associated 
with being in the neigborhood of a quantum critial point 
that can approached via tuning the magnetic-field angle 
with the ah plane llT2l . Reachable within fields H<1Q T, 
the physics of this metamagnetism is acting on a very low 
energy scale of the order of a few meV. Note also that the 
Fermi-liquid regime in vanishing field exists below 10-15 
K, however can be driven to zero temperature with applied 
field 0. 

3 Theoretical Approach In the present modehng 
we combine the LDA approach to DFT with the ro- 
tationally invariant slave-boson theory at saddle-point. 
For the LDA part, we utilized an implementation lfT3l of 
the highly-accurate mixed-basis pseudopotential (MBPP) 
technique |fT4l|. This band-structure code employs norm- 
conserving pseudopotentials jlSl and an efficient mixed 
basis consisting of plane waves and a few localized or- 
bitals. Scalar-relativistic pseudopotentials and the LDA 
exchange-correlation functional after Perdew and Wang llT6l 
were used in the actual computations. 




Figure 1 Left: crystal structure of the orthorhombic Bbcb- 
Sr3Ru207, where the Ru bilayerpart is indicated within the 
dark lines. Right from top to bottom: dyz, d^z and dxy ef- 
fective Wannier function derived from the LDA band struc- 
ture. 



Based on the LDA electronic structure a low-energy 
KS dispersion e^.^^' for each wave vector k within the 
t2g manifold of Sr3Ru207 is constructed via maximally- 
locaUzed Wannier functions (MLWF) 01711181 for the four 
Ru ions in the orthorhombic unit cell. This effective single- 
particle 12 X 12 Hamiltonian is then supplemented by on- 
site Coulomb interactions, resulting in additional Hubbard, 
exchange, spin-flip and pair-hopping terms on each Ru ion 
i with f;2g-orbital indices m, m' in the unit cell a. The in- 
teracting problem hence reads: 

\i.rnm'(T ct 

with the local interacting part in the unit cell given by 
three-orbital generalized Hubbard model on each Ru site, 
i.e., 

im i^m^m' ,(7 

+?7 nira^riim'a + J d\^^d\^^igd^^gd^^i^ 

with n=d^d. The first term in eq. (|2) marks the intra- 
orbital Coulomb interaction with Hubbard U , the second 
term the inter-orbital Hund's rule corrected interaction with 
U'=U—J and U"=U—2J for unequal and equal spin ori- 
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entiation (T=t, i, respectively. Finally, the third term ac- 
counts for spin-flip and the fourth term for pair-hopping 
processes, were we used Jc=J for the real MLWFs. 

For the solution of the problem posed by eq. ([U the 
RISE method was employed, where the electron operator 

dima=dt^a is represented as d^^= RW^X,f^,'a'■ Here R 
is a non-diagonal transformation operator that relates the 
physical operator to the QP operator and is written in 
terms of the slave bosons {0/iri}- In this generalized slave- 
boson theory carries the index A for the physical-electron 
state and n for the QP Fock state. In the end, the general 
idea is to rewrite eq. (|2) solely in terms of slave-boson op- 
erators which allows to integrate out the QP part. The op- 
erator decomposition into QP and interacting slave-boson 
part introduces two constraints, namely 

E'^A„'/'A„ = 1 , (3) 

An 

'l'An'^An{n\flJ,>a^n') ^ fl, f^,,, , (4) 

Ann' 

whereby eq. Q normalizes the total boson weight to unity 
and eq. (HI ensures that QP and boson contents match at 
every site i. This selection of the physical states is imposed 
through a set of Lagrange multipliers {A}. In the mean- 
field version at saddle-point these constraints hold on av- 
erage, with the bosons condensed to c numbers. The nu- 
merical effort amounts to the solution of the saddle-point 
equations for {(/)}, {A} (see Ref. E) for further details). In 
the calculations we also used the set of local Fock states 
{n} for the set of atomic states {A} on each Ru ion. It is 
important to note in this context that in the present work no 
inter-site slave bosons were introduced. Thus at the saddle- 
point, the RISE physical self-energy of the d electrons is 
given by lH 

= (1 - Z-i) + [Rt]-iAR-i - e(o) , (5) 

diagonal in the Ru sites. Here Z=RRT is the QP-weight 
matrix and A is the lagrange-multiplier matrix. The quan- 
tity e'^'^' denotes a possible one-body term in the k-summed 
KS-Hamiltonian, which should not appear in the QP part 
of the framework. Thus Sd contains a term linear in fre- 
quency and a static part. Though more approximative than 
more involved methods like, e.g., quantum Monte-Carlo, 
which may handle the full frequency dependence, this 
proves sufficient in many cases to describe the QP physics 
at low-energy. It also provides additional insight into local 
excitations within a static self-energy approximation. In 
this regard, the slave-boson amplitudes yield direct access 
to the occupation of local multiplets in the metallic state. 
The RISE mean-field approach may therefore be seen as 
meaningful simplified approximation to correlated metals 
beyond LDA. Note that it is superior to LDAh-U in this 
regime, since the latter is designed for Mott-insulating 
systems. 

Our appUcation of LDAh-RISE to Sr3Ru207 considers 
no spin-orbit-coupling (SOC) effects. Albeit SOC is be- 



heved to be important for the low-energy physics of this 
compound, the present study shall exhibit that many inter- 
esting features of the correlated electronic structure may 
already be revealed without that coupling. Further studies 
including the SOC with the obviously well-tailored RISE 
method will be discussed in a future work. 

4 Results and discussion 

4.1 LDA electronic structure The Ru ions in the 
Sr3Ru207 compound are formally in the 4+ state, i.e. the 
id shell of the transition metal is filled with four electrons. 
Since the system is known to be in the low-spin state, all 
these electrons are located in the t2g manifold. Many is- 
sues of the electronic structure on the LDA level have al- 
ready been discussed in previous studies |19 ,20 21J, albeit 
the very low-energy regime has not been addressed theo- 
retically in great detail. 

Figure |2] shows the LDA band structure close to 
the Fermi level, consisting mainly of an Ru(t2g) domi- 
nated 12-band manifold. The dispersion is strongly two- 
dimensional, although still some bands display significant 
variation along F—A in the first Erilluoin zone (EZ). There 
are several rather flat bands close to ep, especially in the 
neighborhood of the X point, as also observed in recent 
angle-resolved photoemission (ARPES) experiments 11211 
l22l . Additionally shown in Fig. [2] are the effective bands 
from the MLWF construction, which are in very good 
agreement with the original LDA bands close to zero en- 
ergy. Due to the strong hybridization of the d^y orbital with 
oxygen states at energies below -2 eV, the effective Wan- 




Figure 2 Erillouin zone with selected high-symmetry 
points (top). The green line marks the chosen path for the 
LDA band structure (bottom) including the downfolded ef- 
fective Wannier t2g low-energy bands (cyan/grey lines). 
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E-Ep (eV) 



Figure 3 LDA density of states normalized to one Ru ion. 
Top: total DOS, middle: local Ru-d DOS within rRu=2.0 
a.u.. Bottom: effective Wannier t2g low-energy DOS. Note 
that the t-zg Wannier functions are tailored to match the 
LDA bands and are thus to some extent build on linear 
combinations of the original d states shown in the middle 
part. 



nier manifold ranges from -2.5 eV to 0.6 eV, i.e., displays 
a bandwidth of about 3.1 eV. 

More details of the LDA spectral behavior of the 
Ru(4(i) states may be extracted from the density-of-states 
(DOS) plot in Fig. [3] It is first observable that the Fermi 
level resides actually between prominent peaks within the 
total DOS of the system. Projecting the local DOS onto 
the standard cubic d harmonics along the cartesian axes 
renders the dominant behavior of the t-zg manifold at low 
energy obvious. Thereby the in-plane d^y orbital has more 
substantial structure close to Ep than the nearly degenerate 
dxz,dyz orbitals. Interestingly, the Cg orbital d^2_y2 has 
also some appreciable weight in that energy regime, with 
important overall clear signatures of hybridization with 
the dxy orbital. The effective Wannier DOS reproduces 
well the low-energy regime, but note that the obtained 
i2g-only Wannier orbitals are now not aligned along the 
cartesian axes anymore, but show some tilting in line with 
the orthorhombic distortions of the Bbcb structure. Fur- 
thermore especially the effective d^y-liks orbital is now 
"dressed" with the dx2-y2 hybridization just mentioned. 
This is indeed visible in the contour plot of the former 
Wannier orbital exhibited in Fig. |4] which also shows the 
slight tilting in the xy plane parallel to the Ru02 bilayer. 
Note that from the underlying MLWF construction the 
on-site energy of the dj;j,-like orbital is higher than for the 
f^az.yz-like by about 1 15 meV. 

4.2 LDA+RISB modeling In order to account for 
electronic on-site correlations beyond LDA we introduced 
the local Hamiltonian Q on each Ru site and solved the 




Figure 4 Isolines of the dxy-H^s Wannier function. The a;- 
and y-axis correspond to the diagonals of the plot. The Ru 
ion associated to the plotted MLWF resides below the (red) 
oxygen ion. 



interacting problem ([T) with the KS t-zg-iiiie Wannier dis- 
persions discussed in the last section. In the following we 
will refer to the latter by d^z.yz.xy, albeit it is understood 
that the effective orbitals are tailored to i35c5-Sr3Ru207 
through the present MLWF construction. The dxz,yz are 
only nearly degenerate, but in the discussion the orbitally 
averaged data is shown, since the present orbital-resolved 
results differ only marginally. Note however that the cal- 
culations allowed for the full differences between all the 
treated local orbitals. According to theoretical estimates 
for the ruthenates II23II24||251 . the value of the Hund's cou- 
pling was set to J=0.35 eV for all the computations. While 
there are LDAh-DMFT investigations for Sr2Ru04 with 
values for the Hubbard U between 1-3 eV II23II26II25L to 
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Figure 5 Top: Orbital-dependent QP weight (top) and or- 
bital occupations per Ru ion (bottom) with increasing U 
for the t2g states. The dotted line in the bottom part marks 
the respective occupation for U =0. 
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Figure 6 Local intra- and inter-orbital spin correlations 
with increasing U . 



our knowledge no such approach exists to Sr3Ru207. Yet 
optics experiments IfZTl point towards a similar U value. 
Here we scanned the onsite Coulomb interaction up to 
U =3.45 eV, since it is also known that due to the neglect of 
quantum fluctuations in slave-boson approaches, the effect 
of U may be underestimated ||28| . 

Figure |5] shows the diagonal QP weight Z„im for 
m=(a;z, yz), xy with respect to U . The inter-orbital terms 
"Zmm' remain marginal and appear irrelevant in the present 
orbital representation. It is seen that moderate U values 
already provide a significant QP renormalization. The lat- 
ter is smaller for the d^y orbital, understandable from the 
larger orbital-resolved bandwidth. In addition. Fig. |5]dis- 
plays the orbital- and JJ-dependent occupations for the to- 
tal four electrons on each Ru ion. As explainable from the 
lower crystal-field level, the dxy orbital is stronger occu- 
pied in the noninteracting case. It may be observed that the 
affect of U on the occupations is rather subtle, with com- 
pensating/polarizing tendencies below/above V* ^235 eV. 
Note of course that the latter value depends strongly on the 
chosen value for J. Though the correlation-induced in- 
terorbital charge transfers are not dramatic, one may still 
expect impact on the very sensitive low-energy physics of 
SraRuaOr. 

According to Hund's rule, in an atomic picture the Ru 
ion should have S= 1 with the dxy orbital being doubly oc- 
cupied on the basis of the given crystal-field splitting. Our 
calculations indeed show that S'^=S(S+V) increases from 
52-1 at J7=0 towards 5^-1.83 for [7=3.45 eV. In this re- 
spect, the orbital-resolved local spin correlations are shown 
in Fig.|6] with the anticipated different behavior for d^y and 
dxz.yz- Namely, S^^ grows strongly with U, whereas 
Sxy diminishes substantially after U* in order to cope with 
the formation of the local spin pair in that orbital. Also the 
inter-orbital spin correlations show the designated xz-yz 
spin-parallel coupling. Another option to investigate the lo- 
cal spin states is given by an inspection of the slave-boson 
weights |0rr'P within the local Ru(t2g) multiplet basis 
{r}. Therefore the slave-boson amplitudes ^„„' are ro- 
tated in the eigenbasis of the isolated local Hamiltonian (|2]i 
via 

<Prr' = l^ln'Pnn'Un'r' > (6) 
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Figure 7 Left: Weights of the local multiplets on a given 
Ru ion via the symmetry-adapted slave-boson amplitudes 
squared with U. Right: Fock decomposition of 

the dominant threefold degenerate triplet state Fj^ |q 
where dxy is the doubly-occupied orbital. 



where lAnr provides the unitary mapping between the mul- 
tiplet basis {r} and the Fock basis {n}. We denote a spe- 
cific multiplet by F^^, where p describes the particle sec- 
tor, r the energy level therein (starting with n=0 for the 
respective ground state) and m marks the spin state, i.e., 
singlet 's', doublet 'd', triplet 't' and quartet 'q'. The diag- 
onal multiplet weights |(/)rrP are plotted in Fig.Qwith re- 
spect to U . Note that in the present case the (j)ppi are close 
to diagonal with only few off-diagonal terms of minor am- 
plitude. It is seen that the local physics is dominated by 
a threefold-degenerate triplet state in the four-particle sec- 
tor with dxy indeed being the doubly-occupied orbital. Be- 
yond that one, two other triplets with only marginal energy 
difference follow and then already a five-particle doublet 
shows up relevant, however loosing weight with increas- 
ing U. 

Besides affecting the local states, the electronic cor- 
relations introduced by U and J of course also have im- 
pact on the QP dispersions. From eq. (|5]l it follows that 
the RISB self-energy narrows the bands via Z and gives 
rise to band shifts through the remaining static term. In 
Fig. [8] we compare the Fermi surface (FS) in LDA (C/=0) 
with the case for [7=3.45 eV. In general the FS is rather 
complex with in total six sheets, labeled 5, ai, 012, 7i, 
and 72 (see e.g. Ref. ||2T]| ). The LDA FS is in rather good 
agreement with ARPES data lISTI . The overall appearance 
bears the prominent quasi-twodimensional character of the 
compound with only weak FS warping along c*, especially 
for the sheets in the inner area of the BZ. However there 
are also bands crossing the Fermi level along c* close to 
the zone boundary, e.g. from M—R. This observation may 
be important for an understanding of the apparent metal- 
lic behavior of the measured optical conductivity along 
c 1291 . The interacting FS shows for the chosen [7 value 
no dramatic differences, which is expected since it is well 
known that LDA provides a surprisingly good fermiology 
for many strongly correlated metals. However there are 
still some important changes that may shed some further 
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Figure 8 Fermi surface of Sr3Ru207 for U=0 (left) and 
C/=3.45 eV (right). In the bottom the view along the c* 
axis in k-space is depicted. 




Figure 9 LDA (dashed lines) and renormalized QP bands 
(J7=3.45 eV) along the high-symmetry directions. The 
lower part shows a blow up around the Fermi level ep with 
the FS sheets 5, ai, a2, 7i, /? and 72. 



light on the generic low-energy physics of Sr3Ru207. The 
renormalized FS can partly also be studied from the devi- 
ations of the associated bands crossing ep along the high- 
symmetry lines in the BZ, as displayed in Fig. |9] For most 
of the sheets there are some minor size changes with U 
(which should be in overall accordance with Luttinger's 
theorem), e.g., 5, ai and 02 somewhat growing in the kz=Q 
plane. However the 72 pocket is rather severely modified. 
First of all from ARPES experiment, 72 appears as a small 
but "simple" hole pocket more or less right inbetween the 
a2 sheet boundary and the X point along the F—X di- 
rection II2TI . On the other hand, the computed (QP) band 



structure in this region looks rather complex, with possi- 
bly as many as four Fermi-level crossings. Intriguingly, 
those crossings pair into two with two bands crossing each 
other in an "dirac-cone-manner", respectively (cf. Fig. |9). 
As it turns out, this low-energy structure is very sensitive 
to electronic correlations and depending on the interaction 
strength, the two named crossing twofolds shift between 
hole- or electron-like fermiology. In the LDA case the 
twofold close to X is higher in energy, giving rise to a hole 
pocket as observed in experiment, whereas the twofold 
close to 02 is barely shifted into the electron-pocket ap- 
pearance. For J7=3.45 eV the situation is nearly reversed 
for the former pocket part, while the latter becomes stabi- 
lized in whats now a sizeable electron pocket. Note that 
since being close to the 72 -pocket structure also the a2 
sheet in the F—X direction is affected by the low-energy 
correlation-induced restructering. This results in a stronger 
renormalization of the a2 sheet in that direction compared 
to the one along F—AF which is in agreement with find- 
ings in recent ARPES experiments Il22l . 

Surely, this whole discussion depends on the magn- 
tiude of U and J, however already the fact that under- 
neath the 72 pocket hides such a correlation-sensitive low- 
energy structure may be of greater importance. The rea- 
son for the sensitivity is besides the flatness of the asso- 
ciated bands close to the X point given by the fact that 
while the two hole-like bands of the outlined fourfold are 
dominantly dxy, the electron-like bands are mainly dxz.yz- 
Thus the multi-orbital character of the compound directly 
manifests in the complex structure around the 72 pocket. In 
the experimental dispersions ETl the described low-energy 
72 structure is partly visible in the occupied area, however 
much more strongly renormalized (lowest-energy band '^l 
meV). The effect of SOC, especially via additional band 
splittings, may here be important to reveal more details of 
this low-energy segment. In addition, the renormalization 
effect through quantum spin fluctuations could be relevant 




-8.15 -0.10 -0.05 0.00 0.05 0.10 0.15 
co(eV) 

Figure 10 RISE QP density of states for the t2g states per 
Ru ion (C/=3.45 eV). Note that the absolute value is some- 
what artificial due to the not-described Hubbard bands in 
the spectral function within slave-boson mean-field and 
should not be directly related to experimental values. 
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in the system since its prone to the ferromagnetic instabil- 
ity. 

Finally, we show in Fig. [TO] the QP DOS from the 
LDA+RISB calculation in a small energy window around 
ep. Albeit we obviously can not fully reach in the present 
modeling the very low-energy scale seen in experiment, it 
is observable that already here a rather subtle energy struc- 
ture shows up in a small range around ep- The latter value 
resides for J7=3.45 eV within a slope of width '-^10 meV 
and onset 6 meV below zero energy. Note that a large num- 
ber of unoccupied states above ep is also verified from 
ARPES measurements ||2T1 . 

5 Conclusions The present work dealt with a real- 
istic LDAh-RJSB approach to the low-energy electronic 
structure of the puzzling Sr3Ru207 compound. It was 
shown that the LDA band structrure close to the Fermi 
level can be accurately downfolded to an effective MLWF 
description with tailored t2g-like orbitals. Supplemented 
by generic on-site Coulomb correlations, the resulting in- 
teracting model accounts for strong electronic self-energy 
effects. A dominant filling of the d^y orbital was identi- 
fied, locally favoring a corresponding four-particle triplet 
state at larger Hubbard U. The renormalization of the low- 
energy electronic structure singles out the 72 pocket (and 
its neighborhood) as especially sensitive to electronic cor- 
relations, also due to its manifest multi-orbital character In 
view of the very low energy scale close to the X point this 
observation draws much attention to that pocket as a possi- 
ble key object for the understanding of metamagnetism in 
Sr3Ru207 Albeit the Hubbard parameters have been 
chosen more or less by hand in this study, the mere fact 
that one may here find definite regions in k-space that are 
rather prominently susceptible to electronic correlations 
seems an interesting result that deserves further detailed 
investigation. The neglect of SOC is a sure drawback of 
the present study, but because of the already good data 
agreement between experiment and theory on, e.g., the 
Fermi surface, very strong changes of the current results 
are not expected with including SOC effects. The latter 
will be of course important when turning on a magnetic 
field H, especially when elaborating on the influence of the 
angle-dependence thereof. However the RISE formalism 
is in principle ideally suited to cope with such physics 
and we plan to address it in a future work. For instance, 
the orbital-occupation shifts in the correlated regime with 
respect to an angle-dependent field, possibly close to the 
metamagnetic response, would be of vital interest. 
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